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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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14.  ABSTRACT 

This  paper  presents  the  extension  of  some  finite  elements  with  embedded  strong  discontinuities  to  the  fully 
transient  range  with  the  focus  on  dynamic  fracture.  Cracks  and  shear  bands  are  modeled  in  this  setting  as 
discontinuities  of  the  displacement  field,  the  so-called  strong  discontinuities,  propagating  through  the 
continuum.  These  discontinuities  are  embedded  into  the  finite  elements  through  the  proper  enhancement  of 
the  discrete  strain  field  of  the  element.  General  elements,  like  displacement  or  assumed  strain  based 
elements,  can  be  considered  in  this  framework,  capturing  sharply  the  kinematics  of  the  discontinuity  for  all 
these  cases.  The  local  character  of  the  enhancement  (local  in  the  sense  of  defined  at  the  element  level, 
independently  for  each  element)  allows  the  static  condensation  of  the  different  local  parameters  considered 
in  the  definition  of  the  displacement  jumps.  All  these  features  lead  to  an  efficient  formulation  for  the 
modeling  of  fracture  in  solids,  very  easily  incorporated  in  an  existing  general  finite  element  code  due  to  its 
modularity.  We  investigate  in  this  paper  the  use  of  this  finite  element  formulation  for  the  special  challenges 
that  the  dynamic  range  leads  to.  Specifically,  we  consider  the  modeling  of  failure  mode  transitions  in 
ductile  materials  and  crack  branching  in  brittle  solids.  To  illustrate  the  performance  of  the  proposed 
formulation,  we  present  a  series  of  numerical  simulations  of  these  cases  with  detailed  comparisons  with 
experimental  and  other  numerical  results  reported  in  the  literature.  We  conclude  that  these  finite  element 
methods  handle  well  these  dynamic  problems,  still  maintaining  the  aforementioned  features  of 
computational  efficiency  and  modularity. 
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The  2007  January  24-26  CMEs 
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First  major  CMEs  observed  by  the  SECCHI  suite... 

...  but  with  a  20  hour  gap  while  the  CMEs  were  in  Mil's  FOV. 

CME\  01/24  @14:03  -700  km/s  -  2^^  CME:  01/25  @06:43  -1200  km/s 
Data  gap:  01/25  06  UT  to  01/26  00  UT 
Studied  by  Harrison  et  al.  (2008)  and  Webb  et  al.  (2009). 

Simulated  by  Lugaz  et  al.  (2009)  and  Odstrcil  et  al.  (2009). 
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STEREO  Behind  C0R1 
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Simulation  set-up 


O  Space  Weather  Modeling  Framework  (Univ.  Michigan) 

O  Solar  wind  model  of  Cohen  et  al.  (2007),  out-of-equilibrium  flux  ropes 
chosen  to  match  initial  observed  speed. 

O  From  IRs,,  to  1  AU:  40,000  4^  blocks  +  15,000  blocks  (>  lOM  cells) 
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Line-of-sight  images 


Simulated  background  is  calculated  by  deriving  the  minimum 
image  from  27  steady-state  LOS  images. 
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Heliospheric  Evolution 


O  Near  solar  nnininnunn:  there  are  a  nunnber  of  steady  structures  (CIR-like). 
O  Interaction  involves  not  only  the  2  CMEs  but  also  these  dense  streams. 


'T=0000  00:00  since  01/24/2007  14UT" 


'T=000000:00  since  01/24/2007  14UT" 


200 


200 


Institute  for  Astronomy 
University  of  Hawaii 


May,  6,2010 


Astronomy 


Heliospheric  Imagers:  what  is 
being  observed? 


2nd  CME  is  brighter  because  the  leading  front  propagates  inside  the  sheath  of  the  1st  CME. 
i^en»e  stream  can  be  identified  as  such  in  these  images,  because  its  propagation  speed  is  different. 
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Time-elongation  plots 


Time  (from  01/25)-elongation  plots  for  PA  69  (apparent  central  PA  of  SECCHI). 
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Simulations  as  laboratory: 
Testing  position  and  velocity  estim^ 
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Here  testing  method  to  derive  position  from  elongation  measurements 
(Kahler  et  al.,  Rouillard  et  al.)- 

id-Phi:  Approximates  CME  with  a  single  point 

^it-P:  Approximates  CME  with  an  expanding  sphere  centered  at  the  Sun 
Position  is  well  determined  up  to  0.5  AU. 

Speed  is  not  very  well  determined  from  existing  methods,  especially  for 
wide,  limb  ejections  past  0.5  AU. 
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New  geometrical  method  to  derive 
CME  position  from  elongation  angle 

(Lugaz  et  al.  Ann  Geo,  2009) 

Instead  of  using  single-point 
approximation  or  assuming  a 
spherically  symmetric  front  we 
use  a  sphere  attached  to  the  Sun 

O  New  assumption  is  good  for  wide 
CMEs  (better  than  Point-P). 

It  can  be  shown  that: 

fp  jj,  =  d  sin  e/  sin  (e-k])) 


r  =  2d  sin  e/(l-t-sin  (E-Kj))) 

l/r  =  .5  (1/ Epp  +1/ Ep 
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Simulations  as  controlled 
laboratory  environment 


O  Testing  fitting  methods  to  derive  CME  direction. 


O  Devised  for  plasma  blobs/CIRs  with  3  assumptions: 

Q  constant  direction,  stereoscopic  methods:  Liu  (2010),  Lugaz  (2010)  in  ApJ 
Q  small  CME  width,  statistical  study  or  different  assumption 
Q  constant  speed.  statistical  study  difficult  (lack  of  fast  CMEs) 


Simulations 

O  Worst  case  scenario  for  fitting  method: 
a  Fast  CME 
Q  Very  wide 

Q  In  front  of  a  coronal  hole  (not  done  here) 


O  We  do  a  simulation  of  a  800  km/s  wide 
CME  (at  20  Rs)  into  a  solar  minimum 
configuration  (bipolar  magnetic  field). 
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Same  CME  from  2  different  view 

points 


O  "Geometrical"  acceleration/deceleration  not  as  clear  as  expected. 

O  Good  for  small  directions  {halo  CMEs),  large  underestimation  of  large 
directions  {limb  CMEs). 

O  MORE  at  my  poster  tomorrow  (EGU2010-3691). 
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°  time 

Best  fit:  606  km/s  and  42° 

[561-694]  and  [27°, 66°] 


°  time 

Best  fit:  619  km/s  and  31° 

[585-659]  and  [22°,45°] 


J-maps  and  direction 


ST-A  75  deg  away 
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Conclusion 


Simulations  are  an  important  tool  in  support  of  HI  observations 
For  complex  events  (CME-CME  interaction): 

Q  What  happened?  (Lugaz  et  al.,  Sol.  Rhys.,  2009,  Webb  et  al,  Sol  Rhys,  09) 
Q  Testing  different  physical  mechanisms  (Lugaz  et  al.,  Ann.  Geo.,  2009) 


O  To  test  and  validate  methods: 

Q  Simulations  are  like  a  controlled  laboratory  environment: 

^CME  speed  and  direction  are  known, 

^  Can  vary  only  one  parameter. 

Q  We  tested  the  Fixed-Phi  and  Point-P  approximations. 

We  devised  a  new  approximation  which  fits  better  our  simulated  CME  and 
1^  is  important  at  large  elongation  angles  (beyond  40°). 

^■Q^ack  of  fast/wide  CMEs  may  explain  why  fitting  method  work  so  well  so 
far. 

O  Thanks  to  NASA  and  NSF  for  funding. 


